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We propose a single-site addressing implementation based on the sub-wavelength localization via 
adiabatic passage (SLAP) technique. We consider a sample of ultracold neutral atoms loaded into 
a two-dimensional optical lattice with one atom per site. Each atom is modeled by a three-level 
A system in interaction with a pump and a Stokes laser pulse. Using a pump field with a node 
in its spatial profile, the atoms at all sites are transferred from one ground state of the system 
to the other via stimulated Raman adiabatic passage, except the one at the position of the node 
that remains in the initial ground state. This technique allows for the preparation, manipulation, 
and detection of atoms with a spatial resolution better than the diffraction limit, which either 
relaxes the requirements on the optical setup used or extends the achievable spatial resolution to 
lattice spacings smaller than accessible to date. In comparison to techniques based on coherent 
population trapping, SLAP gives a higher addressing resolution and has additional advantages such 
as robustness against parameter variations, coherence of the transfer process, and the absence of 
photon induced recoil. Additionally, the advantages of our proposal with respect to adiabatic spin- 
flip techniques are highlighted. Analytic expressions for the achievable addressing resolution and 
efficiency are derived and compared to numerical simulations for 87 Rb atoms in state-of-the-art 
optical lattices. 



I. INTRODUCTION 

Ultracold neutral atoms in an optical lattice with 
single-atom and single-site resolution constitute an ideal 
physical system to investigate strongly correlated quan- 
tum phases [1] which, in turn, has interesting applications 
in quantum optics [2] , quantum simulation [3] and quan- 
tum information processing [4-6], among others. The 
first approaches towards single-site addressing considered 
the use of lattices with relatively large site separations [7] . 
However, to have access to the regime of strongly corre- 
lated systems, typical lattice spacings well below 1 /jm are 
needed since the tunneling rate has to be comparable to 
the on site interactions. In this case, the diffraction limit 
imposes strong restrictions on the addressability of indi- 
vidual lattice sites. To overcome this limitation different 
techniques have been investigated. For instance, spatially 
dependent electric and magnetic fields have been used to 
induce position dependent energy shifts on the atom [8], 
allowing for site-selective addressability. Alternatively, 
a scanning electron microscopy system to remove atoms 
from individual sites with a focused electron beam [9] has 
been reported. However, in this case, atoms need to be 
reloaded into the emptied sites after each detection event. 
More recently, high resolution fluorescence imaging tech- 
niques, that make use of an optical system with high nu- 
merical aperture, have been implemented to perform in 
situ single-atom and single-site imaging for strongly cor- 
related systems [10]. In this context, a single-site address- 
ing (SSA) scheme based on focused laser beams inducing 
position-dependent energy shifts of hypcrfine states has 
been theoretically [11] and experimentally reported [12]. 
In the experiment, an intense addressing beam is tightly 



focused by means of a high resolution optical system. 
This beam produces spatial dependent light shifts bring- 
ing the addressed atom into resonance with a chirped 
microwave pulse and eventually inducing a spin-flip be- 
tween two different hyperfine levels of the atom. 

On the other hand, during past years, several propos- 
als based on the interaction of spatially dependent fields, 
e.g., standing waves, with three-level atoms in a A-type 
configuration have been considered, not only for single- 
site addressing in optical lattices but, more generally, for 
sub- wavelength resolution and localization [13-16]. In 
the first approaches [13-15], a spatially modulated dark 
state is created by means of either electromagnetically 
induced transparency (EIT) or coherent population trap- 
ping (CPT) [17, 18], which allows for a tight localization 
of the atomic population in one of the ground states, 
around the position of the nodes of the spatially depen- 
dent field. More recently, it has been shown that the 
resolution achieved with those CPT or EIT based tech- 
niques can be surpassed using stimulated Raman adia- 
batic passage (STIRAP) [19] processes, by means of the 
so-called sub-wavelength localization via adiabatic pas- 
sage (SLAP) technique [16]. The SLAP technique relies 
on a position-dependent STIRAP of atoms between the 
two ground states of a three-level A atomic system, and 
has additional advantages compared with CPT or EIT 
techniques such as (i) robustness against parameter vari- 
ations, (ii) coherence of the transfer process, that allows 
for its implementation also in Bose-Einstein condensates 
[16], and (iii) the absence of photon induced recoil. 

In this work, we apply the SLAP technique [16] to ul- 
tracold atoms in an optical lattice, where single-site ad- 
dressing (SSA) requires one to overcome the diffraction 
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limit. In order to address only a single site, we use here 
Stokes and pump pulses with Gaussian shaped spatial 
distributions, with the pump presenting a node centered 
at the lattice site that we want to address. Assuming 
that all the atoms in the optical lattice are initially in 
the same internal ground state and applying the stan- 
dard STIRAP counterintuitive temporal sequence for the 
light pulses [19], we will demonstrate that it is possible to 
adiabatically transfer all the atoms, except the one at the 
node of the pump field, to an auxiliary ground state. We 
will show that this process is performed with higher effi- 
ciency and yields better spatial resolution than the CPT 
based techniques [13-15]. Also, we will demonstrate that 
our addressing technique requires shorter times than in 
the adiabatic spin-flip technique discussed in Ref. [12], 
and that larger addressing resolutions can be achieved 
using similar focusing of the addressing fields. Moreover, 
our technique has the additional advantage that it can 
be applied between two degenerated ground-state levels. 

The paper is organized as follows. In Sec. II we in- 
troduce the physical system under consideration. In Sec. 
Ill, we present a protocol to achieve SSA, and we derive 
analytical expressions for the spatial resolution and ad- 
dressing efficiency of our technique. In addition, a com- 
parison with CPT based techniques is provided. Next, 
in Sec. IV, we perform a numerical investigation of the 
proposed technique for a single-occupancy optical lattice 
loaded with 87 Rb atoms by integrating the corresponding 
atomic density-matrix equations. Finally, in Sec. V, we 
summarize the results and present the conclusions. 



II. PHYSICAL MODEL 

Figure 1(a) illustrates the physical system under con- 
sideration. It consists of a sample of ultracold neutral 
atoms loaded into a two-dimensional (2D) square optical 
lattice with spatial period A/2, placed in the plane (x, y) 
and illuminated by a pump and a Stokes laser pulse with 
Rabi frequencies fip and Og, respectively, propagating 
in the —z direction with a selectable time delay. The 
spatial profile of the pump pulse has a node coinciding 
with a particular lattice site, our target site, at which 
the Stokes pulse is also centered. With the spatial pro- 
files of the pulses having revolution symmetry around the 
propagation axis [dashed line in Fig. 1(a)], in the follow- 
ing we consider only the transverse spatial dimension x, 
without loss of generality. In our model, we assume the 
system to be in the Mott insulator regime with only one 
atom per lattice site. Each atom is considered to have 
only three relevant energy levels in a A-type configura- 
tion, defined by the interaction with the light pulses as 
shown in Fig. 1(b). Here, 721 (723) is the spontaneous 
transition rate from the excited state |2) to the ground 
state |1) (|3)) and A P (Ag) is the detuning of the pump 
(Stokes) field. We assume that all atoms are initially in 
state |1). 

Our approach to achieve single site addressing is based 




FIG. 1. (Color online) (a) Physical system under investiga- 
tion: the pump and Stokes light pulses, with Rabi frequencies 
Qp and Qs, propagate in the —z direction and interact with 
the atoms of a single-occupancy optical lattice located in the 
(x, y) plane, (b) Scheme of the A-type three-level atoms, ini- 
tially in state |1), that interact with pump and Stokes pulses. 
Excited level |2) has spontaneous transition rate 721 (723) 
to level |1) (| 3}) and Ap (As) is the detuning of the pump 
(Stokes) field. 

on the SLAP technique [16], where, depending on their 
position, the atoms are transferred between two internal 
ground states by means of the STIRAP technique [19]. 
STIRAP consists in adiabatically following one of the 
energy eigenstates of the A system, the so-called dark 
state, which under the two-photon resonance condition, 
i.e., Ap = Ag, has the form 

\D(x,t)) = cos8(x,t) |1) - sin0(x,i) |3) , (1) 

where tan8(x,t) = Qp(x,t)/Qs(x,t). Starting with all 
the population in |1), it is possible to coherently transfer 
the atomic population to state |3) changing adiabatically 
8 from 0° to 90° by means of a convenient time sequence 
of the fields. This time sequence corresponds to apply 
first the Stokes pulse and, with a certain temporal over- 
lap, the pump pulse. Since the process involves one of 
the eigenstates of the system, the population transfer is 
robust under fluctuations of the parameter values if these 
are adiabatically changed and the system does not evolve 
near degenerate energy eigenvalues. 

In the SLAP technique, the pump field has a spatial 
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structure with nodes yielding state-selective localization 
at those positions where the adiabatic passage process 
does not occur, i.e., those atoms placed at the nodes 
of the pump field remain in 1 1) , while those interacting 
with both fields, pump and Stokes, are transferred to 
|3). For our purposes, we use the SLAP technique with 
a pump field having a single node at the position of the 
target site. Therefore, at the end of the SLAP process the 
population of all atoms illuminated is transferred from 1 1) 
to 1 3) except for the one at the node of the pump field. 
The spatial and temporal profiles for pump and Stokes 
Rabi frequencies are given by 
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where flpo and flso are the peak Rabi frequencies, tp 
and ts are the centers of the temporal Gaussian profiles, 
wp and ws are the spatial widths of the node in the 
pump and of the Stokes field, respectively, and a is the 
temporal width. 

There exist several methods to create the required 
pump intensity profile with a central node: e.g. (i) re- 
imaging of a Gaussian beam with a dark central spot cre- 
ated by a circular absorption mask, using (ii) a Laguerre- 
Gaussian laser beam [20] or (iii) a "bottle beam" created 
by the interferometric overlap of two Gaussian beams 
with differing waists [21], or (iv) a flexible intensity pat- 
tern generated by spatial light modulators and subse- 
quent imaging [22, 23]. 



III. SINGLE-SITE ADDRESSING 

In our model, we assume that the spatial wavefunctions 
of the individual atoms placed at the different sites, cen- 
tered at x n (being n the site index), correspond to the 
ground state of the trapping potential, which in first ap- 
proximation can be considered harmonic. Therefore, the 
full atomic distribution in the lattice is given, initially, 
by 
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where u> a t = \Jh/muj is the width of the initial atomic 
distribution at an individual site, m is the mass of the 
trapped atom and w is the harmonic trapping frequency. 
We assume that the addressed site is xq = and their 
nearest neighbors x±i are at a distance ±A/2 where A is 
the wavelength of the fields that create the optical lattice. 

In order to characterize our single-site addressing tech- 
nique we consider that, once the SLAP technique has 
been applied, the final atomic population distribution in 



|1), pf LAP (x), is given by 

,SLAP/ 



pr—'W = Pgt P (*)Pux(x), (5) 

where P^ AP (x) is the probability distribution that an 
atom remains in state |1) after the SLAP process. Using 



the SLAP technique, the addressing resolution that one 
can obtain is related to the global adiabaticity condition 
[19] at each spatial position x, 



soe 



-x 2 /wj 



PO 



-X j Wj 



> 



(6) 



where T = tp — ts and A is a dimensionless constant 
that, for optimal Gaussian temporal profiles and overlap- 
ping times, takes values around 10 [19, 24]. In Eq. (6), 
the equality gives a spatial threshold Xth above which 
the adiabaticity condition is fulfilled. Assuming that the 
full width at half maximum (FWHM) of Pj^f(x) is 
(Ax)slap ~ %th and expanding Eq. (6) up to first order 
in x one obtains 
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where R 1 = 
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tion (7) gives the width of the addressing region, and it 
tends to zero as R increases. Moreover, since (Ax) SLAP 
must be real valued, we find that the inequality 
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must be fulfilled. In this paper we will consider that Q so 
is fixed, so R' can be varied through fipo, wp, and ws- 

Two conditions should be satisfied for our SSA tech- 
nique to work. First, the population of the atom in the 
addressed site must remain in state |1) after the action 
of the fields and, second, the rest of the atoms of the lat- 
tice have to be transferred to level |3). Therefore, taking 
into account the overlap between P-^ AP (x) and pi a t(x) 
it is clear that the FWHM of the probability 
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distribution P^f p (x) should satisfy 



(Ax) a t < (A£)slap < x x - (Ax) at , 
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where (Ax) at = 2\/ln 2 iy at , and xi — A/2 is the position 
of the nearest-neighboring site. Using Eq. (7), it is easy 
to see that these conditions fix the range for flsoT to 
obtain SSA using the SLAP technique: 
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with x+ = (Ax) a t and x_ = x\ — (Ax) a t. Note that the 
upper limit for Eq. (10) is more restrictive than Eq. (8). 

In order to have a quantitative description of the SSA 
performance, let us introduce the SSA efficiency as 



v = Vx a (l-V Xl ), 



(12) 
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where V Xo corresponds to the probability of finding the 
atom at the addressed site Xq in state |1), while 1 — V Xl 
corresponds to the probability that the atom in the neigh- 
bor site X\ has been transferred to a different internal 
state. We define 
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with ±s = Xi ± A/4 and i = 0,1, whereas pf LAP (x) and 
piat (a;) have been defined in Eqs. (5) and (4), respectively. 
Using Eqs. (4) and (5), the explicit forms for Eqs. (13) 
are 
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From these expressions, it can be seen that, for x\ 
(Aa;) a t, the limits given by Eq. (9), i.e., (Ax) SLAP = 
(Aa;) at and (Ax) SLAp = X\ — (Ax) at , correspond to SSA 
efficiencies of ^0.70 and ^0.94, respectively. 

An alternative technique to perform atomic localiza- 
tion based on spatial dependent dark states is the co- 
herent population trapping (CPT) technique [13]. In the 
CPT [17] technique the dark state is populated after sev- 
eral cycles of coherent excitation followed by spontaneous 
emission from |2) to the ground states. Note that, while 
CPT relies on spontaneous emission, the SLAP technique 
is fully coherent. Moreover, the latter provides higher 
resolution, as shown in Ref. [16], and does not suffer from 
recoil since the localized atoms have not interacted with 
light. In what follows, we compare the range of param- 
eters necessary to perform SSA considering both SLAP 
and CPT techniques. Note that we focus our compar- 
ative analysis in the CPT technique, although similar 
results are obtained considering the method proposed 
in Ref. [15], where the spatial dependent dark state, 
\D(x, t)), is created via the STIRAP technique by switch- 
ing off the fields before completing the transfer process. 

In order to compare both techniques, we define the fi- 
nal population distribution in 1 1) using 

CPTaspfPT^) 

in an analogous way as it has been done for the SLAP 
technique in Eq. (5). Then the FWHM of the correspond- 
ing probability function P^^(x) is obtained by impos- 
ing that |(1 \D(x,t))\ 2 = 1/2 and t P = t s in Eqs. (2) and 
(3): 
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Fixing the desired Ax and using Eqs. (7) and (16) for 
SLAP and CPT, respectively, the constraints for the rel- 
evant parameters for each technique can be obtained. For 
simplicity, we consider that the Stokes pulse parameters 
ws and flso are fixed, and only the node width wp and 
pump peak Rabi frequency flpo can be varied. Note that, 
for the SLAP case, we have to fix also A and T. 
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FIG. 2. (Color online) (a) Parameter values of R needed to 
perform SSA as a function of wp, using SLAP (solid line) and 
CPT (dashed line) techniques. The addressing probability 
distribution widths are taken as Ax — A/4 in both cases, (b) 
(Aa;)sLAP (solid line) and (Ax)cpt (dashed line) as a function 
of R. The FWHM of the atomic distribution, [Ax) at = A/10 
corresponds to the horizontal dotted line and we have taken 
nip = Ai. The parameters used in both (a) and (b) are 
ws = 32Al = 24A, Q. S0 T = 19, and A = 20. 



Taking Ax = A/4, half of the site separation, the re- 
quired values for R and wp are plotted in Fig. 2(a). To 
simultaneously illuminate a large number of sites, we use 
a large Stokes beam waist of ws — 32 Al = 24 A, where 
\l = | A is the wavelength of both pump and Stokes 
fields [25]. The solid (dashed) line corresponds to the 
SLAP (CPT) case with the parameter values QsoT = 19 
and A = 20 [24]. As icp decreases, in both SLAP and 
CPT cases, lower values of R are needed to reach the 
fixed resolution Ax, since the narrower the node of the 
pump field, the narrower the probability distribution of 
atoms remaining in |1). In addition, for any given width 
of the node of the pump field w p , the required values of 
R are lower in the SLAP case than in the CPT case. 

It is important to realize that Eqs. (7) and (16) 
show the possibility to obtain, for certain parameter val- 
ues, widths of the probability distribution, (Ax)slap or 
(Ax)cpt, smaller than (Ai) at . In particular, using the 
SLAP technique this can be achieved with moderate R 
values. This is shown in Fig. 2(b), where (Ax)slap (solid 
line) and (Ax)cpt (dashed line) are represented as a 
function of R for wp = Al and the rest of the parameters 
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as in Fig. 2(a). The FWHM of the atomic distribution, 
(Ax) a t = A/10, is depicted with a horizontal dotted line 
to indicate the values where (Ax)slap < (Ax) at . As we 
stated in the discussion of Eqs. (14) and (15), this limit 
corresponds to a SSA efficiency of 77 ~ 0.70. This regime 
of parameters is interesting because it shows that the 
SLAP technique could be used for applications in site- 
selective imaging with a resolution down to the width of 
the atomic distribution at each site. 



IV. NUMERICAL SIMULATIONS 

In this section, by numerically integrating the corre- 
sponding density-matrix equations, we study the imple- 
mentation of the SLAP-based SSA technique for A-type 
three-level 87 Rb atoms in a single-occupancy optical lat- 
tice. Numerical calculations using the CPT technique are 
also presented for comparison. The wavelength of the 
lasers that create the optical lattice, red detuned with 
respect to the Dl line of 87 Rb, is A = 1064 nm. The po- 
tential depth of the optical lattice is chosen as Vb = 15E r , 
where E r is the recoil energy. This corresponds to a har- 
monic trapping frequency of w = 2-7T x 15.92 kHz [3]. 
Therefore, the FWHM of the atom distribution at each 
site due to the confining potential is (Ax) a t = 142 nm. 
Pump and Stokes fields with A^, = 795 nm are cou- 
pled to |1) -H> 1 2) and |3) f-> |2), respectively, where 
|1) = |F = 2,m F = -2), 1 2) = \F' = 2,m F = -1) and 
1 3) = \F = 2, mp = 0) are hyperfine energy levels of the 
Dl line of 87 Rb. The excited state |2) has a spontaneous 
transition rate 721 = 2n x 0.96 MHz (723 = 2n x 1.44 
MHz) to state |1) (|3)), and we assume no spin deco- 
herence during the interaction time. We consider equal 
temporal pulse widths of a = 0.2 /is with a temporal 
delay T = 1.4c, in such a way that the total SSA pro- 
cess time is 4cr. The Stokes pulse has a maximum Rabi 
frequency fiso = 19/T = 2n x 10.8 MHz, while the max- 
imum Rabi frequency of the pump is varied through the 
parameter R, since fipo = flsoVR- Concerning the spa- 
tial profiles of the fields, we assume a wide Stokes profile, 
ws = 32 A^, and a narrow node for the pump, wp = Aj> 
As it has been discussed in the previous section, to prop- 
erly perform SSA, the population of all the atoms in the 
lattice, except the one in the addressed site, must be 
transferred from |1) to |3) with high probability. In what 
follows, those requirements for the realization of the SSA 
are studied by numerically evaluating the SSA efficiency. 

The signatures of SSA are shown in Fig. 3, where the 
numerically evaluated efficiency and probabilities defined 
in Eqs. (12) and (13) are plotted as a function of R for 
both the SLAP (circles) and CPT (crosses) techniques. 
Figures 3(a) and 3(b) show the probabilities V XQ and 
1 — V Xl , respectively. For large values of R the prob- 
ability of finding the atom at the addressed site in |1) 
is higher with the CPT than with the SLAP technique 
[see Fig. 3(a)], while for small values of R the probabil- 
ity of removing the atom from the neighboring site [see 



1 

0.8 
o 0.6 
k 0.4 

0.2 



0.2 


1 

0.8 
0.6 





6 


1 "x 




X 




X 




X X X X 

° o 


x X x XXXX) 


- 















- 




... 1 






1 


(a) " 


cr 
- 





'"6 




X 




X 


6 6 ' '6 

v x x 
x x 


° § 999m 
- 






X 










- 


X 










- 


- 












(b) " 


cr 















— x— x— xaooas 






,-'''x 




G 




-° 




X 
























(c) " 



1 10 100 1000 

R 



FIG. 3. (Color online) Numerical results for the probability of 
finding the atom located at xo in state |1) (a), the probability 
to transfer it from |1) to another state (b), and the efficiency 
») as a function of R (c), for the SSA with SLAP (circles) and 
CPT (crosses) techniques. Analytical curves for SLAP (solid 
line) and CPT (dashed lines), computed from Eqs. (12) and 
(13), are added in (c) for comparison (see text for the rest of 
the parameters). 

Fig. 3(b)] is higher in the SLAP case. This is explained 
by the fact that (Ax) slap < (Ax)cpt for a given value of 
R, as shown in Fig. 2(b). The SSA efficiency results ob- 
tained by the numerical evaluation of Eqs. (12) and (13) 
are shown in Fig. 3(c), together with the corresponding 
analytical curves (solid and dashed lines) obtained us- 
ing Eqs. (14) and (15), added for comparison. A good 
agreement is found between numerical and analytical re- 
sults. From Fig. 3(c) it is clear that the SLAP technique 
is more efficient for lower values of the intensity ratio of 
the addressing fields (R < 50) than the CPT technique. 
Certainly, this is advantageous for the experimental im- 
plementation with limited laser power available. 

An example of the final population distribution after 
performing SSA with the SLAP (circles) and the CPT 
(crosses) techniques is plotted in Fig. 4(b) for the par- 
ticular case of R = 10 and with the rest of the param- 
eters as in Fig. 3. The initial population distribution 
in |1), pi a t(x), at the addressed site (xo = 0), and the 
two next neighbors x±\ — ±532 nm is shown as a solid 
line. Note that, in the SLAP case (circles), the popu- 
lation of state |1) around xq remains almost the same 
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FIG. 4. (Color online) (a) Spatial profile of the pump pulse 
with its node focused at x = 0. (b) Final population dis- 
tribution remaining in |1) using SLAP (circles) and CPT 
(crosses) single-site addressing techniques for R = 10. The 
initial atomic distribution in pi at (x), is shown as a solid 
line (see text for the parameters values). 



after the addressing process, while in the first neighbor 
sites it is practically zero. On the other hand, for the 
CPT case (crosses), the population in the addressed site 
remains also nearly unchanged, but it exhibits a signifi- 
cant amount of population in the neighbor sites. This is 
in full agreement with the discussion following Figs. 3(a) 
and 3(b). For this example, the total efficiencies found 
are 77 CPT = 0.56 and 7y SLAP = 0.95 according to the corre- 
sponding values shown in Fig. 3(c). In addition, as shown 
in Fig. 4(a), the width of the pump node required to per- 
form the SSA method is much larger than the addressed 
region (Ax) slap- In particular, for the case shown in 
Fig. 4, wp = \l and (Ax)slap — ^l/3, thus obtaining 
addressing resolution beyond the diffraction limit. 

Finally, we carry out a comparison between our pro- 
posal and the experiment reported in Ref. [12], where 
a focused laser beam induces position-dependent light 
shifts, allowing one to perform a spin flip by means of 
a resonant microwave pulse at the addressed site. Since 
the microwave field involved in the experiment has a Rabi 
frequency of kHz, the total spin- flip time is in the order 
of ms. In contrast, as our proposal makes use of only op- 
tical fields, the addressing time is three orders of magni- 
tude below (ps) . Specifically, to achieve similar values of 
the addressing resolution in both techniques, (Ax) ~ 300 



nm, we have obtained an addressing time of ~40 fis. This 
decrease of the total addressing time needed implies a re- 
duction of the effects caused by spontaneous scattering 
of photons, which are a limitation for the light shifts- 
based proposals [11, 12]. These effects could be strongly 
reduced in our case. Also, we have compared the res- 
olution obtained in both techniques. In Ref. [12], they 
use an addressing beam with an intensity FWHM of ap- 
proximately 600 nm, and obtain a spin-flip probability 
distribution with FWHM=330 nm. In our technique, us- 
ing R — 1 and a width of the pump node wp = 509 
nm, which corresponds to the width of their addressing 
beam, a very similar FWHM of the addressing proba- 
bility distribution is obtained: (Acc)slap = 330.66 nm. 
Note that this value can be reduced by increasing the 
ratio between the pump and the Stokes intensities, e.g., 
R = 10 implies (A:e)slap = 181.86 nm and for R = 100 
we obtain (Ax)slap = 100.82 nm. 



V. CONCLUSIONS 

In this work we have discussed a proposal to perform 
single-site addressing (SSA) in an optical lattice by us- 
ing the SLAP technique. With respect to other dark 
state based techniques such as CPT, this method is fully 
coherent, robust against variations of the parameter val- 
ues, and we have found that it yields higher efficiencies 
for smaller values of the intensity ratio between the pump 
and the Stokes fields. Moreover, the addressed atom does 
not interact with the fields, minimizing all possible de- 
coherent processes. On the other hand, with respect to 
the recent experiment on SSA using adiabatic spin flips 
[12], the present proposal allows one to use two degen- 
erate ground levels, takes shorter times to perform the 
addressing process, and provides similar or even larger 
addressing resolutions for similar focusing of the address- 
ing fields. 

The proposed method provides an achievable address- 
ing resolution that can be pushed well below the diffrac- 
tion limit of the addressing light field and of the optical 
setup used for addressing or detection of atoms at closely 
spaced lattice sites. This relaxes the requirements on the 
optical setup or extends the achievable spatial resolu- 
tion to lattice spacings smaller than accessible to date. 
Through analytical considerations, we have derived the 
range of parameters for which SSA is properly achieved. 
Moreover, we have obtained expressions to estimate the 
resolution and the efficiency of the SLAP-based address- 
ing method, and we have compared them with the anal- 
ogous expressions obtained using the CPT technique. 
Next, by integrating the density-matrix equations with 
realistic parameter values for state-of-the-art optical lat- 
tices loaded with 87 Rb atoms, we have checked the valid- 
ity of the analytical approach. 
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